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D ■ Abstract 

Homothetic scalar field collapse is considered in this article. By making 

in _ 

C*~> . a suitable choice of variables the equations are reduced to an autonomous 

O 

o> 

O 

\ shown that there are two classes of solutions. The first consists of solutions 

q-i with a non-singular origin in which the scalar field collapses and disperses 

(3JT), again. There is a singularity at one point of these solutions, however it is not 



X 



system. Then using a combination of numerical and analytic techniques it is 



visible to observers at finite radius. The second class of solutions includes 

both black holes and naked singularities with a critical evolution (which is 

neither) interpolating between these two extremes. The properties of these 

solutions are discussed in detail. The paper also contains some speculation 

about the significance of self-similarity in recent numerical studies. 
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I. INTRODUCTION 



That gravitational collapse leads to black hole formation is widely accepted, yet com- 
paritively little is known about the generic features of collapse. This is well exemplified 
by the lack of a precise formulation of the cosmic censorship hypothesis ffl]. The reason 
for this is simply our lack of the mathematical tools necessary to analyse the evolution of 
generic initial data, although some progress is being made on this front Moreover the 
complexity of the Einstein field equations counsels retreat to a more tractable system which 
may capture the essence of the physics. Spherically symmetric general relativity coupled to 
a massless scalar field is one such system. Indeed this model has been studied in great detail 
both analytically |§ and numerically [£|||. Christodoulou has rigorously established global 
existence and uniqueness of solutions to the Einstein-scalar field equations and has discussed 
the general properties of these solutions 0. He has even established a sufficient condition 
for the formation of a trapped surface in the future evolution of a given initial data set ||. 
Numerical investigations have also provided useful insights into black hole formation. The 
most recent of these studies, by Choptuik f7|, has revealed several intriguing phenomena 
which were hitherto unknown. 

Choptuik considered the numerical evolution of initial data sets characterised by a single 
parameter (p say). The resulting families of solutions, S[p], include both geometries contain- 
ing black holes, and geometries with only slight deviations from flatness, depending on the 
value of p. Between these two extremes lies a critical evolution with p = p*, which signals the 
transition between complete dispersal and black hole formation. His most interesting results 
pertain to near critical, p ~ p*, evolutions which exhibit a particularly simple strong field 
behaviour. In fact, his results strongly suggest that near critical evolutions may be described 
by a single, universal solution of the field equations. Two quantitative features also emerge 
from his work: (i) Near critical evolutions contain echoes in the strong field region (i.e. form 
invariant quantities exhibit a scaling relation like ( |1.2j ) below, with A ~ —3 • 4). Based on 
his extensive investigations Choptuik has conjectured that the exactly critical solution has 
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an infinite train of echoes in the strong field region approaching the singularity, (ii) When 
black holes form in near critical evolutions the mass scales as ~ \p — p*\@, where (3 is 
an apparently universal constant determined numerically to be about -37 

Obtaining analytical results about solutions which exhibit a discrete self-similarity of the 
type discussed by Choptuik has proved to be extremely difficult [IIJ. In ]TTJ it has been 
argued that the piling up of the echoes as the singularity is approached suggests that a 
continuous self-similarity may be a good approximation in this neighbourhood. Therefore it 
seems worthwhile to examine collapse under the assumption of a continuous self-similarity. 
Furthermore I would like to draw attention to a possible explanation of the echoing, dis- 
covered by Choptuik, in terms of a continuous self-similarity. In curvature coordinates the 
spherical line element may be written as 

ds 2 = -a 2 dt 2 + a 2 dr 2 + r 2 dtt 2 (1.1) 

where a and a are functions of both r and t. Self-similarity implies that there exists a time 
coordinate such that the metric depends only on the combination x = r/t, consequently 
quantities like a(r, t) satisfy the relation 

a(e A r,e A t) = a(r,t) (1.2) 

for arbitrary values of A. Now, consider a coordinate tranformation to a new time, T say, 
in terms of which t = f(T) where the function / satisfies 

f{e A T) = e A f{T) (1.3) 

only for discrete values of A. In terms of this new time coordinate the form invariant 
quantities examined by Choptuik, along with a(r/ f(T)), would all appear to have a discrete 
(rather than continuous) self similarity. How could this be the origin of the discreteness in 
the numerical work, after all Choptuik J?] states quite clearly that he uses central proper 
time to describe the solutions - this corresponds to the similarity time t in a spacetime 
which is everywhere self-similar. However it is likely, and indeed suggested by numerical 
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work, that (discrete) self-similarity holds only on intermediate scales, < t\ < r < r2 say. 
(This has been further highlighted in numerical studies of radiation fluids Thus it is 
possible that the time coordinate used in the numerical studies does not correspond to the 
similarity time. It should be noted that even if this is true the appearance of the universal 
scaling constant A ~ 3 • 4 continues to be a mystery. 

In any case, this paper is concerned with homothetic spacetimes; it is assumed that there 
exists a vector field £ such that the spacetime metric satisfies £%g = 2g, where denotes 
the Lie derivative with respect to £. The homothetic collapse of perfect fluids has received 
a great deal of attention during the past few years |I2| , where the main thrust of work 



was searches for naked singularities fL3| . The problem of homothetic, scalar field collapse 
has received comparitively little attention ||, and no examples of naked singularities which 
evolved from regular initial data were previously known for this form of matter. 

The reduction of the scalar field equations to an autonomous system allows me to give a 
detailed description of all the solutions. Two interesting features emerge from the analysis: 
Firstly the existence of solutions with naked singularities, and secondly the occurence of 
phase transitions in the system. The critical point behaviour is of two distinct forms. The 
first is a phase transition from solutions which (roughly) represent dispersal, to geometries 
representing black holes. It has been suggested elsewhere [14| that the critical solution 



implied by the numerical results may be both on the verge of black hole formation and 
being a naked singularity. We will see below that the second phase transition corresponds 
to such a situation - the critical evolution lies at the boundary between black holes and 
naked singularities. In this article the term black hole is used rather loosely to mean that 
an apparent horizon exists and precedes a central singularity. Of course it is possible to 
obtain an asymptotically flat solution to the Einstein-scalar equations by cutting off the 
self-similar evolution at some advanced time and matching it to a less symmetric (not self- 
similar) exterior, in this way one would obtain an asymptotically flat, black-hole spacetime. 

The paper is arranged as follows: In section II the field equations for the self-similar 
collapse are derived in terms of a retarded time coordinate u and a radial coordinate r. 



The most general (see appendix A) scalar field evolution consistent with the homothetic 
symmetry of the spacetime is ip — h(r/\u\) — k\d.\u\, where k is an arbitrary (positive) 
constant. These equations have already been derived in |J and used to provide initial data 
in a search for naked singularities. A field redefinition transforms the equations into a non- 
linear autonomous system which is amenable to standard techniques ||16|| . It is then shown 
that the solutions fall into two distinct classes depending on the value of k. Solutions of the 
first type (k>1/ \^4tc) are shown to be singular only at one point on r = 0. This singularity 
is not naked; distant observers cannot see it by waiting a finite proper time. The transition 
point is at k — 1/ V^tt, although the exactly critical solution actually belongs to the first 
class. 

For each k < I/a/Itt there exists a continuous infinity of solutions with a non-singular 
origin. For a given value of k solutions can represent black holes or naked singularities, with 
a critical point evolution interpolating between these two extremes. The analytic examples 



of phase transitions which have been discussed in the literature |TT| are of this type, in fact 
they correspond to the case k = when the equations can be integrated in closed form. 
This solution is briefly discussed to explain how it fits into the general analysis presented 
here. Section IV contains some discussion of the results, and their possible significance. 



The notation of |L5f is adopted throughout the paper, and detailed calculations which 
might distract from the main line of thought are relegated to the appendices. 



II. THE FIELD EQUATIONS 

Using retarded Bondi coordinates {u, r, 8, 0} the spherical line element may be written 

ds 2 = -ggdu 2 - 2gdudr + r 2 dtt 2 (2.1) 

where g = g(u } r),g — ~g{u, r) and dTt 2 is the standard line element on the unit two-sphere. 
The origin is singular unless ~g = g when r = 0. Furthermore it is convenient to normalise 
the coordinate u so that it represents proper time for an observer at the origin, thus write 
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g(u,0) = g(u,0) = l. (2.2) 

The Einstein-scalar field equations are then 

(ln^) ir . = 47rr(^) 2 , (2.3) 

(rg),r = 9 , (2.4) 

(g/g),u = [(ip, u ) 2 ~ gip,uip,r] , (2.5) 

where a comma denotes partial differentiation, if) = ip(u, r) is a massless, minimally coupled 
scalar field satisfying 

(gr 2 ^,r), r = 2rV, u + 2rV,r« . (2.6) 
Spherical symmetry allows the introduction of a local mass function m(x a ) jl9| defined 



by 



1 - := 9 »"r„r s = I , (2.7) 





where r is the function which determines the area of the two-spheres. This mass function 
agrees with both the ADM and Bondi masses in the appropriate limits, and is equivalent to 
the Hawking quasi-local mass in this case. 

A. Self-similar ansatz 

The existence of a homothetic symmetry in a spherical spacetime implies that the metric 
depend only on x = r/\u\, and that the scalar field evolve as 

^ = h(x) -Kki\u\ (2.8) 

where h is some function to be determined and k is constant (see Appendix A for a proof 
of this fact). Writing 

h {x) = r lip. di ( 2.9) 

Jo 4 
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the self-similar equations derived from (|2.3|) - ([2.5|) are 



W = 9 , 

xg' = ^ng{if , 
g - g = ^ \2k 2 x - (g - 2x)(^ 2 + 2Kj) 



(2.10) 
(2.11) 
(2.12) 



where a prime (') denotes differentiation with respect to x. In deriving (|2.12| ) from (|2.5|) it 
is necessary to use ( |2.10| ) and ( |2.11| ) to eliminate derivatives of g and g. The scalar field 
evolution is determined by 



x(g — 2x)j' = 2kx — / y(g — 2x) . 



(2.13) 



It is now a straightforward matter to show that ( 2.10 )-( 2.12j) imply (|2.13|) provided 7 7^ —K. 
These equations have been derived previously by Goldwirth and Piran 0] and used to provide 
boundary conditions in a numerical search for naked singularities. 

At the origin ( j2.2|) and ( 2.12j ) imply either 7(0) = or 7(0) = —2k. Directly evaluating 
the trace of the stress-energy tensor for the scalar field one finds 
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7 2 - =(7 2 + «7) 
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,0 v 



(2.14) 



Clearly the solution can have a non-singular origin only if 



7(0) = . 



(2.15) 



This completes the specification of the initial conditions for the equations (|2.10|) - ([2.12|) . It 
is however convenient to recast the system in an autonomous form before discussing the 
solutions. 



B. An equivalent autonomous system 



Analysis of the above equations is facilitated by the field redefinitions 

y = g/g, z = x/g, 
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(2.16) 



and the introduction of a new coordinate 

f = lnx. (2.17) 



Upon substitution into ( 2.10| ), ( |2.11| ) and ( 2.13 ) one obtains the three dimensional, non-linear 
autonomous system 



z = z[2-y- 1 ) , (2.18) 
ij = 1 - (4tt 7 2 + l)y , (2.19) 
(l-2z)7 = 2/?z- 7 (?r 1 -2z) . (2.20) 

The system is effectively two dimensional, however, since 7 is determined by the algebraic 
relation 



7 = —K ± 



^ 4tt(1-2z) ' 

provided 2/7^ 1/(1 + Attk, 2 ). Further discussion is therefore couched in terms of a projection 
into the yz-plane. Consistent with the initial condition ( |2.15| ) the positive square root is 
taken in ( |2.21[ ), however it must be emphasised that solutions may still evolve continuously 
onto the other leaf of the surface defined by taking the negative square root above. Indeed 
it is solutions of this type which have naked singularities. 

Requiring that the mass function, defined in (|2.7| ), should be positive or zero implies 
y < 1, while 7 is real only if 

(1 + A-IXK 2 ) - y- 1 

± - — — > . 2.22 

l-2z K 1 

Black hole formation is signalled by y — > (technically this is the condition which locates 
an apparent horizon in the spacetime). The continuation of the solution to negative values 
of y will not be considered in the sequel, thus, integral curves of interest lie in the strip 
< y < 1. 

It should be noted that 7 is not continuous at 

,2 = 1/2, y = 1/(1 + 4ttk 2 ) , (2.23) 



thus invalidating the usual existence and uniqueness theorems for systems of ordinary dif- 
ferential equations at this point. This has the important consequence that integral curves 
of the differential equations may intersect at Q2.23| ). Furthermore, the continuation of such 
solutions is not always uniquely defined, in some cases there exists an infinite family of 
possibilities. This is discussed in more detail below and in appendix C. 

The character of the solutions shows a strong dependence on the value of k. Solving 
( p.lO| )-( [2~T3| ), subject to the regularity conditions on r = 0, gives rise to two essentially 
different classes of solution according as 4ttk 2 is greater than or less than unity; they are 
solutions which do not contain black holes, and solutions which contain either black holes 
or naked singularities. The critical evolution (4ttk, 2 = 1) also belongs to the first class. The 
transition from one class to the other as one adjusts k is similar to the behaviour discussed 
by Choptuik 0. A second type of phase transition from solutions containing black holes 
to those with naked singularities occurs for each value of k in the range < 47r/t 2 < 1. 
This may simply be an artifact of the restriction to self-similar continuations past the point 
( |2.23j ) however no other continuations are considered here. A similar phenomenon has been 
observed in Tolman-Bondi collapse flT^|. 

The system ( |2.18| )-( |2T2Ti| ) has two stationary points, one on either leaf of ( |2.21| ), given by 



The nature of these points depends on the value of k. It is discussed below and in Ap- 
pendix B. 



III. SELF-SIMILAR SOLUTIONS 

The self-similar solutions fall quite naturally into two distinct classes depending on the 
value of k. Class I solutions do not contain black holes or naked singularities although they 
are singular at one point on r = 0. In the second class are solutions which have naked 
singularities or black holes; a single critical evolution having a null singularity (which is not 
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naked) interpolates between these two extremes. For completeness, at the end of this section 
it is shown how the results of [II] fit into the overall picture. 



A. Class I - Attk 2 > 1 

For k in this range only the stationary point (y + , z+, 7+) is of interest. This exact solution 
is equivalently written as 



g = 2x(l + VAttk) , g = x(l + VAttk) . (3.1) 

It has a singular origin (r = 0), in fact there are two sheets of this singularity. Future 
directed, ingoing lightrays terminate on the sheet located at u — 0, while outgoing lightrays 
originate on the past sheet (The solution may be obtained by setting a = (3 = in Eq. (9) 
of 0). 

Locally, (|3TT|) is a positive attractor (see appendix B). The global structure of the Class I 
solutions is easily determined by examining the behaviour of the integral curves in the region 

A = {1/(1 + 4vr/t 2 ) <y<l, 0<z< 1/2} . (3.2) 

Noting that integral curves enter A across the lines y — 1/(1 +4ttk, 2 ), y = 1 and {z = 0, y > 
1/2} it is evident that the solution originating at z — 0, y — 1 either terminates at the 
stationary point or leaves A across z — 1/2. In appendix C it is shown that integral curves 
only cross z = 1/2 at y = 1/(1 + 47tk 2 ), and that the solution passing through this point is 
unique when 47tk 2 > 1. A direct consequence of this is that the solution with a non-singular 
origin approaches the asymptotic form (|3.1| ), its evolution being characterised by a sequence 
of decaying oscillations in y about the value y = 1/2. (See Fig. 1) 

These solutions do not contain black holes, yet there is a singularity at the origin when 
u — > 0. Is this singularity naked? Consider ingoing, radial null geodesies 

dr , . . 

- = -3/2 . (3.3) 

Since x = —r/u the geodesies are given by 
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y 

FIG. 1. Integral curves when 4ttk 2 = 2.25. The solutions are seen to spiral towards the 
stationary solution. The dot-dashed curve from y = 1, z = is the solution with a non-singular 
origin. The dashed integral curve is on the negative leaf of 7 and is the continuation of the solid 
curve which reaches y = 1/(1 + Attk 2 ). 
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\n\u\ = 2 / (g-2x')- l dx' , (3.4) 

where Xi is the radius at the point of intersection between the ingoing geodesic and the null 
cone u = —1. Since z = x/~g < 1/2 throughout the entire evolution the integrand in ( |3.4| ) 
is always bounded, and x decreases with increasing u. Thus ingoing lightrays must reach 
r = before u — 0, and they never reach the singularity provided X\ < oo. 

In view of the asymptotic solution ( |3.1| ) one also finds that an observer at a (large) fixed 
radius takes an infinite proper time to reach u = 0. A central observer, on the other hand, 
reaches the singularity at u = in finite proper time. It is, however, only after he has seen 
the entire history of the universe in a tremendous flash. Thus these spacetimes have trivial 
topology, with a singularity only at r = as u — > 0. 

Similar arguments apply when Atck 2 = 1, although in this case ~g — 2x — > as i ^ oo. 
It must be emphasised that this solution is unstable, in the sense that an arbitrarily small 
change in the value of k drastically changes the character of the resulting solution. In 
particular, for smaller values the spacetime contains a black hole. 



B. Class II - 4vtk 2 < 1 

Containing two sub-classes of solutions and offering another example of a phase transition 
in gravitational collapse - from black hole spacetimes to naked singularities - these solutions 
are more interesting than those in Class I. In fact for each k in the range < Atck 2 < 1 
there appears to be a continuous infinity of solutions with a regular origin. This is due to 
the failure of uniqueness at the singular point (|2.23|) . Naked singularities develop only for 
sufficiently small values of k. 

Solutions with a non-singular origin contain a null hypersurface, T say, on which x = con- 
stant. This corresponds to the point y = l/(l+47Tfv 2 ), z = 1/2 in phase space. As mentioned 
earlier the standard uniqueness theorems break down at this point, and there is a one param- 
eter family of self-similar continuations beyond it. More general extensions which produce 
asymptotically flat spacetimes have been considered in p. 
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FIG. 2. Integral curves when 4ttk 2 = 0-49. Here we see the integral curve representing the 
solution with a non-singular origin reach the point z = 1/2, y = 1/(1 + Attk 2 ). There is a one 
parameter family of continuations past this point. All solutions of interest terminate at y = (an 
apparent horizon in spacetime) except one which approaches the stationary point. Here also the 
dashed lines are integral curves on the negative leaf of 7. 
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FIG. 3. A spacetime diagram for the critical evolution when Attk 2 < 1. The unshaded region 
corresponds to the dot-dashed curve in Fig. 2 and Fig. 4. The shaded region shows a continuation 
past T [where z = 1/2, y = 1/(1 +4ttk 2 )] which approaches a stationary point. The singularity at 
r = is null. 
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Fig. 2 shows various solutions to the system ( |2.18| )-( p.20| ) when 4ttk 2 = • 49. It is clear 
from the diagram that all solutions (except one) evolve into black holes (signalled by the 
formation of an apparent horizon, y — > 0). Some of the solutions evolve onto the negative 
leaf of ( |2.21| ); y undergoes a single oscillation before decreasing monotonically to y — 0. 
While I have no analytic proof that all of these solutions contain black holes, the numerical 
results (as illustrated in Fig. 2) suggest this is true. The single exceptional solution exhibits 
a behaviour which has been discussed elsewhere ; beyond T no black hole forms, instead 



the solution asymptotically approaches one of the stationary points (|2.24|) . The spacetime 
is singular at the null surface u = 0, r = where y = 1/2. This singularity lies at infinite 
redshift for observers at large radius. Fig. 3 is a spacetime diagram for this exceptional case. 

For smaller values of k more complicated behaviour is possible, and indeed another type 
of phase transition is apparent - from black holes to naked singularities. When z > 1/2 
the x = constant surfaces are spacelike. If no black hole forms then x — > oo as u — > 
[see Eq. ( |3.4j )1, corresponding to a Cauchy horizon in the spacetime. For sufficiently small 
k there are solutions which evolve from the singular point onto the negative leaf of ( ^.21[ ), 



asymptoting to y = 1/(1 + Attk, 2 ) as z — > oo. A typical solution of this type is shown in 
Fig. 4, for Attk 2 = • 25. It is not difficult to obtain an approximate solution in the large 
z limit, and hence to show that the Cauchy horizon is non-singular. For large z Eq. Q2.19|) 
becomes 

y ~ 1 - (1 + AixK 2 )y . (3.5) 



Integrating and inserting the result in ( 2.181 ) implies 



^TTirf +CT "" +W) (3 ' 6) 

z = % - ~ dx 1 -*™ 2 , (3.7) 
9 

where c and d are arbitrary constants of integration. It is now straightforward to verify the 
regularity of the Ricci scalar on the Cauchy horizon, using fl2.21|) and ( 2.14Q 



K 2 
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FIG. 4. Integral curves when 4ttk 2 = 0-25. 77 is an example of an integral curves which evolves 
onto the negative leaf of 7 but does not reach y = 0, instead it attains a minimum value of y and 
then asymptotes to y = 1/(1 + 4itk 2 ). Such solutions have a naked singularity at r = 0. 
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Only as r — > is this quantity singular, indicating the existence of a naked singularity at 
u = 0, r = in these solutions. The Cauchy horizon is also a null orbit of the homothetic 
Killing vector, therefore the work of Lake and Zannias implies that the singularity is 
strong in the sense that tidal forces diverge at it. Notice also that the metric can be made 
manifestly regular at the Cauchy horizon by transforming to the new coordinate given by 



dU = (-w)- 4 ™ du . (3.9) 

These spacetimes confirm that naked singularities, which evolve from regular initial data, 
also exist for scalar field sources. Whether they are stable to non-homothethic, not to 
mention non-spherical, perturbations is an open question although the work of Goldwirth 
and Piran || suggests that they are not. 

There are two critical points where phase transitions are observed. Each of the exactly 
critical solutions, which interpolates between naked singularities and black holes, asymptot- 
ically approaches one of the stationary points ( |2.24 ). Therefore these spacetimes have the 



structure shown in Fig. 3, with a null singularity at u = 0, r = 0. 

As n decreases further, the topology of the phase space changes slightly; only a sin- 
gle phase transition occurs going from black holes to naked singularities. The solutions 
(including the single critical evolution) are the same as those already discussed. 

C. The Roberts solution k = 



When k = the scalar field is a function of x = —r/u only and the equations (|2.10|) - (|2.12|) 



can be exactly integrated [fnj| . The resulting solution was first discovered by Roberts f20 
and may be written as 

g = {(p 2 + x 2 ) 1/2 -2p 2 }/x (3.10) 



g = 2g[l±y/l- Ap 2 (2xg - g 2 ) ) (3.11) 

and 



17 



1.0 



0.8- 



0.6- 



0.4- 



0.2- 



0.0 




0.0 0.2 0.4 0.6 0.8 1.0 

y 



FIG. 5. The exact expression for the solutions when k = is given in the text, however this 
diagram is for comparison with the other cases. Note that y=l is now a solution to the equations 
and corresponds to flat space. For this reason no naked singularities exist for any value of p when 
k = 0. 
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7 = p/^4:7t(p 2 + x 2 ) . (3.12) 

The integral curves are plotted in Fig. 5. The solutions are labelled by p and exhibit critical 
point behaviour. When p — 1/2 the integral curve leaves z — 1/2, y — 1 and approaches the 
stationary point fl2.24| ). This exactly critical evolution lies between solutions which contain 



black holes (y — > 0, z — > 1/2 along the integral curves) and those which evolve back to flat 
space (y — > 1/(1 + Ann 2 ), z — > oo along the integral curves). Since this is a saddle point we 
see that near critical evolutions (sub- or super-critical) can approach this point arbitrarily 
closely before moving away in their respective directions. 

k = is therefore exceptional because the subcritical evolutions have zero mass on the 
Cauchy horizon, u = 0, as is readily seen by taking the limit x — > oo in the above solution. 
This behaviour is also apparent in the asymptotic approach of the integral curves to the line 



y = 1. As discussed in |TT| the natural continuation past u = is Minkowski space since 
there is no material flux across this surface on which m = 0. Notice however that there is 
no self-similar extension of these solutions to r = which is non-singular. 



This discussion places the results obtained in [O]] within the more general context of 
self-similar spacetimes with scalar field matter sources. 



IV. DISCUSSION 

Spherically symmetric, homothetic spacetimes have received a great deal of attention 
over the last few years due to the ease with which it is possible to construct naked singu- 
larities in such spacetimes. Recent numerical studies of spherical collapse suggest that 
self-similarity may play an important role in describing the approach to the singularity in 
gravitational collapse. This study of scalar field collapse was in fact motivated by the work 
of Choptuik, where he observed discrete self-similarity in solutions on the verge of black hole 
formation. However it is not clear how a continuous self similarity, as discussed here, could 
be at the center of the results which he has obtained. Nevertheless some interesting features 
do emerge from the study of spacetimes with a homothetic symmetry. 
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The scalar field evolves as ip = h(r/\u\) — K,ln\u\ where h(r/\u\) is determined by the 
coupled Einstein- scalar field equations. The constant k distinguishes between two different 
classes of solution. The first class (when 47tk 2 > 1) are non-singular almost everywhere. 
Technically there is a singularity at u = = r, however it is not naked and the worldlines 
which reach it in finite proper time are a set of measure zero. When < Attk 2 < 1 the 
solutions are more interesting - for each value of k in this range there exists infinitely 
many solutions with a non-singular origin. When 4nK 2 is only slightly less than unity all 
of these spacetimes contain apparent horizons. However for sufficiently small n some of the 
spacetimes have naked singularities. 

47tk 2 = 1 marks the transition point from solutions of Class I to those of Class II - in 
effect from no black holes to black holes. Another phase transition occurs for each value of 
k < 1/ V^r - a transition from black holes to naked singularities. In particular the critical 
evolution is on the verge of being a naked singularity in this case. One might wonder if there 
is some way to classify the critical point behaviour more generally, and to which such class 
the generic case belongs. 

One of the most intriguing results obtained by Choptuik was the scaling law for black 
hole mass, unfortunately self-similar spacetimes cannot have finite mass black holes. In order 
to obtain an asymptotically flat spacetime it is necessary to cut off the self-similar evolution 
at some advanced time, and consider a suitable continuation (which is not self-similar). 
Goldwirth and Piran || did exactly this, although they did not examine the behaviour 
of black hole mass as the critical point 47m 2 was approached. This question is currently 



under active investigation [[1711 , and it will be interesting to see if the mass exhibits the same 
behaviour which has been observed elsewhere 

We have in self-similar scalar field collapse further examples of spacetimes which violate 
cosmic censorship. One might be surprised about this were it not for the plethora of examples 
which now exist. What emerges from these examples (generally) is that there do exist initial 
data sets which when evolved according to the Einstein equations lead to naked singularities, 
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however the genericity of these data is far from clear. It therefore seems that the thrust of 
any attempt to formulate (and prove) cosmic censorship must address this issue directly. 



Some interesting preliminary results have been obtained by Lake J^j where he has shown 
that (spherically symmetric) spacetime in the neighbourhood of a naked singularity may be 
approximately self-similar. It therefore seems that future work on naked singularities must 
consider deviations from the symmetric situations treated to date. 

Finally in searching for a theoretical understanding of the results obtained by Chop- 
tuik and Abrahams and Evans || one might consider the obvious generalisation of home- 
theticity to a conformal symmetry; that is to suppose the existence of a vector field £ such 
that 

C 6 g = n(x)g . (4.1) 

If the dependence on position in Q is weak, solutions might behave like self-similar solutions 
with some sort of super-imposed periodicity. It would seem interesting to investigate this 
possibility. 
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APPENDIX A: SCALAR FIELD EVOLUTION IN SPHERICAL, SELF-SIMILAR 

SPACETIMES 

Results of Defrise-Carter |23[] imply that a spherical spacetime with a homothetic sym- 
metry (i.e. there exists a vector £ such that C^g = 2g) can be written in the form 

ds 2 = e 2t ( gi (x)dt 2 + g 2 (x)dx 2 + e 2x dVl 2 ) (Al) 
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where dVt 2 = dO 2 + sin 2 9 dip 2 , and the similarity vector is £ = d/dt. 

The exponential dependence of the metric on t gaurantees that the Christoffel symbols 
and hence the Ricci tensor, R^, are independent of this coordinate. This may be expressed 
covariantly as 

CtRp, = . (A2) 

Before examining the implication of this for a self-similar spacetime which satisfies Einstein's 
equations with scalar field matter, let me show that (|A1|) can be recast into the form used 
in section II. 

Introduce new coordinates r and u defined by 

r = exp(t + x), u — rG(x), (A3) 

where G(x) is to be determined. Substituting them into the line element ( |AI| ) and requiring 
u to be null one obtains the ordinary differential equation 

d ^ = G(-l±^-g 2 {x)/g l {x)) (A4) 

which determines G(x). Furthermore the line element reduces to 

ds 2 = —g{r/u)g(r/u)du 2 — 2g(r / u)dudr + r 2 dVl 2 , (A5) 

where r/u is realted to x by ( |A3|) . Clearly this means that g and g are functions of r/u as 
stated in section II. The similarity vector is 

£ = rd r + ud u (A6) 

in these coordinates. 

Now in view of (|A2| ) and the Einstein field equations 

R P v = 8vn/vVv ( A7 ) 

the scalar field must satisfy 
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= • (A8) 



(A9) 
(A10) 



Assuming that ip is independent of 6 and 0, Eq. ( |A8[ ) is readily integrated to 

dip 7(r/n) 

<9r r 

<9-?/> — j(r/u) 

du u 

where 7 is an arbitrary function of r/u. (Integrability was used to reduce the number of 
arbitrary functions to one) The general solution of these coupled equations is therefore 

ip — h(r/u) — kIu \u\ — /31n |r| . (All) 



In particular (3 may be set to zero by absorbing (3(ln \u\ — In |r|) into h, thus reducing ( All|) 

to am. 



APPENDIX B: THE STATIONARY POINTS 



The discussion in section III relies on the properties of the stationary points of the 



equations (|2.18 )-( 2.20 ). In general there exists two such points given by 



v ± = i, z± 



1 



1 



7± 



(Bl) 



1 ± y/InK ' \ ' \~ 

Provided 47tk 2 7^ 1 it is straightforward to linearise about each of these points and hence 
to analyse the topology of the phase space in their neighbourhoods. The eigenvalues are 



Ai, 



2z 4 



± 



l - z ± f + 4(2z ± - i; 



(B2) 



1 \ (2^ ± -l) 2 

where z± may be chosen independently of the sign of the square root. 

When Attk, 2 > 1 only z + is relevant for the discussion in section III. Both eigenvalues are 
real and have the same sign, Ai,2 < 0, when 1 < 4tik 2 < 4/3 so that the stationary point is 
an attractive node. Once 47tk 2 > 4/3 the eigenvalues become complex conjugate, and since 



2z + - 1< , 1 - z+ > 



(B3) 
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they have negative real part. This is a positive attractor, with spiral behaviour. 
When < Attk, 2 < 1 both stationary points are of interest. Simply noting that 

2z ± - 1 > , 1 - z+ > , 1 - 2_ < (B4) 

the discussion of both is easily combined. Clearly the eigenvalues are real and have opposite 
signs since 



|1 - z ± \ < y/(l - z ± ) 2 + A{2z ± - 1) . (B5) 
Thus they are saddle points. 

APPENDIX C: THE SINGULAR LINE Z = 1/2 

In the above analysis it is important that integral curves cannot cross z — 1/2 except at 
y = ot at y = 1/(1+ Attk 2 ). We now show that this is so, and derive the solution in the 
neighbourhood of z = 1/2, y — 1/(1 + Attk 2 ). The analysis is split into two cases: 

(i) Suppose an integral curve crosses z = 1/2 at y ^ 1/(1 + Ann 2 ). Writing z = 1/2 + ( 
and considering the ( — > limit of ( |2.18| ), ( |2.19| ) and ( [2.21| ) it is a straightforward matter to 
derive 

dy ^ (1 + 4TTK 2 )y 2 - y 

d(~ (2y-l)C ' 1 1 

Integrating this equation gives 

(1 — A-KK 2 ) 

In (CI ~ constant + t rf In 1(1 + AtckV - y\ . (C2) 

11 (1 + 47TK 2 ) ' V ' ' K ' 

Examining this expression shows that y^ = is the only place where integral curves may 
intersect z = 1/2. 

(ii) It is necessary to treat the case when CI J/ =i/(i+4 7r/c 2) = separately since the limit z — > 1/2 
in (|2.21|) is more delicate. For this purpose we introduce ( as above and write 
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where 77 <C 1/(1 + 47tk 2 ), thus 



i + 4ttk 2 

7~-k± A — '- . (C4) 



-57T 



-c 

Substituting this approximate expression for 7 into (|2.19|) , and using ( as the independent 
variable we arrive at the equation 



77V 2 ~ ± - — -L - 9A (C5) 



rfC 2C(1-47tk 2 )' (-C) 1 / 2 (1-4ttk 2 ) 



and find 



,,1/2 ~ ^(-C) 1/2 + c (- C )( 1 + 4 - 2 )/(^ 4 - 2 ) , (C6) 

where c is a constant of integration. Now the initial condition is 77 = when ( = so that 
there are two distinct possibilities; when 4ttk, 2 > 1 the integration constant must vanish, 
implying that a single integral curve passes through this point. On the other hand if 47tk 2 < 1 
the constant is not fixed by the initial conditions and there is a one parameter family of 
curves passing through 77 = = £■ It is exactly this fact which gives rise to the variety of 
solutions in Class II. 
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